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I. INTRODUCTION 



An analytical solution of the radial part of the Klein-Gordon equation is of high im- 
portance spinless in relativistic quantum mechanics, because the wave function contains all 
necessary information for full description of a quantum system. There are only few poten- 
tials for which the radial part of the Klein-Gordon equation can be solved explicitly for all 
n and I. So far, many methods were developed, such as supersymmetry (SUSY) [1,2] and 
the Pekeris approximation [3-8], to solve radial part of the Klein-Gordon equation exactly 
or quasi-exactly for I ^ within these potentials. 

The one- dimensional Klein-Gordon (KG) equation is investigated for the PT -symmetric 
generalized Woods-Saxon (WS) potential [9] and Hulten [10] and is solved by using the 
Nikiforov-Uvarov (NU) method which is based on solving the second-order linear differential 
equations by reduction to a generalized equation of hypergeometric type. 

The radial part of the Klein-Gordon equation for the Woods-Saxon potential [11] can not 
be solved exactly for 1^0. It is well known that the Woods-Saxon potential is one of the 
important short-range potentials in physics. And this potential were applied to numerous 
problems, in nuclear and particle physics, atomic physics, condensed matter and chemical 
physics. Therefore, it would be interesting and important to solve the radial of the Klein - 
Gordon equation for Woods-Saxon potential for 1^0, since it has been extensively used to 
describe the bound and continuum states of the interacting systems. 

In this work, we solve the radial part of the Klein-Gordon equation for the standard 
Woods-Saxon potential using NU method [12], and obtain the energy eigenvalues and cor- 
responding eigenfunctions for any I states. 

II. NIKIFOROV-UVAROV METHOD 

The Nikiforov-Uvarov (NU) method is based on the solutions of general second-order 
linear equations with special orthogonal functions. It has been extensively used to solve 
the non-relativistic Schrodinger equation and other Schrodinger-like equations. The one- 
dimensional Schrodinger equation or similar second-order differential equations can be writ- 
ten using NU method in the following form: 



where a(z) and a(z) are polynomials, at most second-degree, and t(z) is a first-degree 
polynomial. 

Using Eq(2.1) the transformation 

4>(z) = *(z)y(z) (2.2) 
one reduces it to the hypergeometric-type equation 

a(z)y" + r(z)y' + Xy = 0. (2.3) 



The function <&(z) is defined as the logarithmic derivative [10] 

&(z) n(z) 



(2.4) 



$(z) <t(z)' 
where ■k(z) is at most the first-degree polynomial. 

The another part of tp(z), namely y(z), is the hypergeometric-type function, that for fixed 

n is given by the Rodriguez relation: 

where B n is the normalization constant and the weight function p(z) must satisfy the con- 
dition [12] 

±(a(z)p(z))=T(z)p(z), (2.6) 

with t(z) = t(z) + 2n(z). 

For accomplishment of the conditions imposed on function p(z) the classical orthogonal 
polynomials, it is necessary, that polynomial t(z) becomes equal to zero in some point of an 
interval (a, b) and derivative of this polynomial for this interval at a(z) > will be negative, 
i.e. t'(z) < 0. 

The function 7r(z) and the parameter A required for this method are defined as follows: 



.2 

o — T . fa' — T 



«( z ) = ——±^——j -a + ka, (2.7) 

\ = k + n'(z). (2.8) 

On the other hand, in order to find the value of k, the expression under the square root 

must be the square of a polynomial. This is possible only if its discriminant is zero. Thus, 
the new eigenvalue equation for the Eq.(2.1) is [12]: 

A = \ n = -nr> - I^-ll a", (n = 0, 1, 2, ...), (2.9) 
After the comparison of Eq.(2.8) with Eq.(2.9), we obtain the energy eigenvalues. 



III. SOLUTIONS OF THE KLEIN-GORDON EQUATION WITH THE WOODS- 
SAXON POTENTIAL 



The standard Woods-Saxon potential [11] is defined by 

V(r) = - V ° (r - RoV a«flo, 0<r<oo. (3.1) 
1 + exp (__a) 

This potential was used for description of interaction of a neutron with a heavy nucleus. 
The parameter i? is interpreted as radius of a nucleus, the parameter a characterizes thick- 
ness of the superficial layer inside, which the potential falls from value V = outside of a 
nucleus up to value V = — V inside a nucleus. At a = 0, one gets the simple potential well 
with jump of potential on the surface of a nucleus. 

In the spherical coordinates, the stationary Klein-Gordon equation with Woods-Saxon 
potential is [13] 



1 r 2 dr 1 dr ' 



r 2 sin 6 m \ Sm dl) + r 2 sin 2 9 d^ 2 



i(j=[(E- V{r)f - m 2 c 4 ] if). 

(3.2) 



where m is the rest mass of a scalar particle. 

The terms in the square brackets with the overall minus sign are the dimensionless angular 
momentum squared operator, L 2 . Defining ip(r,9,{p) = R(r)Y(9,<p), we obtain the radial 
part of the Klein-Gordon equation [13] 



d 2 R{r) 2 dR{r) ^ 



dr 2 



r dr 



(e - vy - my + i 



h 2 c 2 



R{r) = 0. 



(3.3) 



where I is the angular momentum quantum number. 
After introducing the new function 



u(r) = rR(r), 



Eq.(3.3) takes the form 



d 2 u(r) 
dr 2 



(E - VY - m 2 c 4 1(1 + 1 



u(r) = 0. 



(3.4) 



h 2 c 2 r 2 

The Eq.(3.4) has the same form as the equation for a particle in one dimension, except 
for two important differences. First, there is a repulsive effective potential proportional to 



the eigenvalue of h 2 l(l + l). Second, the radial function must satisfy the boundary conditions 
u(0) = and u(oo) = 0. 

If in Eq.(3.4) introduce the notations 

r- Ro R 
R a 

then the Woods-Saxon potential is given by the expression 

V - V ° 

VWS — -T- ■ -, r- 

1 + exp(o;xJ 

It is known that the radial part Klein-Gordon equation cannot be solved exactly for this 
potential at the value I ^ using the standard methods as SUSY and NU. From Eq.(3.6) 
it is seen, that the centrifugal potential V/(r) is inverse square potentials, which cannot be 
solved analytically. Therefore, in order to solve this problem we can take the most widely 
used and convenient for our purposes Pekeris approximation. This approximation is based 
on the expansion of the centrifugal barrier in a series of exponentials depending on the 
internuclear distance, taking into account terms up to second order, so that the effective I 
dependent potential preserves the original form. It should be pointed out, however, that this 
approximation is valid only for low vibrational energy states. By changing the coordinates 
x = ^—j^ 2 - or r = Rq(1 + x), the centrifugal potential is expanded in the Taylor series around 
the point x = (r = R ) 

= w + 1) = w + 1) i _ _ 

lK J 2m r 2 2m R 2 (1 + x) 2 V ; ' V ' 

where 5 = % . 

According to the Pekeris approximation, we shall replace potential V/(r) with expression 
[7,8]: 



vr (r) = 5 ( C + -r—^: + - ; 2 - 2 , (3.6) 



C\ i C<2 

1 + expax (1 + expax)' 
where Co, C\ and C2 are some constants. 

In order to define the constants Co, C\ and C2, we also expand this potential in the Taylor 
series around the point x = (r = Rq): 



vr(x) = s 



„ C\ C 2 \ a a 2 2 a 3 ~ a 4 ^ 4 

Co + + ^ - T (A + £2) £ + tttCsx 2 + — (Ci + C 2 ) x 3 - — C 2 x 4 + 



2 4 / 4 v y 16 48 v y 96 

(3.7) 



Comparing equal powers of x Eqs.(3.5) and (3.7), we obtain the constants Co, C\ and C2 
[7,8]: 

4 12 8 48 48 
^0 — 1 1 — w — ^2 — 9 • 



a cr 



cr 



Instead of solving the radial part Klein-Gordon equation for the centrifugal potential 
Vi(r) given by Eq.(3.5), we now solve the radial part Klein-Gordon equation for the new 
centrifugal potential VJ*(r) given by Eq.(3.6) obtained using the Pekeris approximation. 
Having inserted this new centrifugal potential into Eq.(3.4), we obtain 



dr 2 



+ 



, , (2EV «i+i)a \ ( vj 1(1+1)0* 

E 2 -mQC A l(l + l)C \ \h 2 c 2 a 2 a 2 J \h 2 c 2 a 2 a 2 



h 2 c 2 



a 2 a 2 



1 , 

1 + e " 



(l + e^)' 



u = 0. (3.8) 



We use the following dimensionless notations 



e = — 



(E 2 - m 2 c A ) a 2 I (I + 1) C \ 2 _ 2EV a 2 I (I + 1) C x _ 2 



h 2 c 2 



cr 



h 2 c 2 



cr 



; 7 = - 



V 2 a 2 Z (/ + 1) C 2 



/t, 2 c 2 



(3.9) 



CT 



with real e > (E 2 < m^c 2 ) for bound states; j3 is real, 7 is real and positive. 
If we rewrite Eq.(3.8) by using a new variable of the form 

-1 



1 + exp 



we obtain 



... , 1 — 2z ,. . — e 2 + B 2 z — 7 2 z 2 . . 
u"(z) + — -u'(z) + — f = 0, (0 < 2 < 1), 



z(l-z) v 7 (z(l-z)) 2 

with f (z) = 1 - 2z; a(z) = z{\ - z); a(z) = -e 2 + fi 2 z - 7V. 
In the NU- method the new function n(z) is 

n{z) = ±^e 2 + (k - (5 2 ) z - {k - 7 2 ) z 2 . 



(3.10) 



(3.11) 



The constant parameter k can be found employing the condition that the expression 
under the square root has a double zero, i.e., its discriminant is equal to zero. So, there are 
two possible functions for each k 

( e - Je 2 -f3 2 + 7 2 ) z-e, for k = f3 2 - 2e 2 + 2eJe 2 - (3 2 + 7 2 , 

n(z) = ±{ > V / V (3.12) 

(e + y/e 2 - /3 2 + 7 2 J z - e, for k = (3 2 - 2e 2 - 2e^/e 2 - /3 2 + 7 2 . 
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According to the NU-method, from the four possible forms of the polynomial ■k(z) we 
select the one for which the function t(z) has the negative derivative and root lies in the 
interval (0,1). Therefore, the appropriate functions tv(z) and t(z) have the following forms 

n(z) = t -[e+^e 2 -p 2 + 7 2 ) z, (3.13) 



t(z) = 1 + 2e - 2 (l + e + ^e 2 - /3 2 + 7 2 ) z, (3.14) 

and 



k = P 2 - 2e 2 - 2ey/e 2 - /3 2 + 7 2 . (3.15) 
Then, the constant A = k + 7r'(;z) is written as 



A = (3 2 - 2e 2 - e 2 - /3 2 + 7 2 - e - ^e 2 - (3 2 + 7 2 . (3.16) 
An alternative definition of X n (Eq.(2.9)) is 

A = \ n = 2 (e + ^e 2 - /? 2 + 7 2 ) n + n(n + 1). (3.17) 
Having compared Eqs.(3.16) and (3.17) 

/3 2 -2e 2 -2e v / e 2 - /3 2 + 7 2 -e- ^ t 2 - (3 2 + 7 2 = 2 (e + ^e 2 - /3 2 + 7 2 ) ra+ra(ra+l), (3.18) 
we obtain 



/ i \/i + 4 7 2 

e + x/e 2 - /? 2 + 7 2 + n + - - V . r = 0, 



(3.19) 



2 2 
or 

e+ Ve 2 - /3 2 + 7 2 - n' = 0. (3.20) 

Here 

„,_ n + yTW-i i (321) 

n being the radial quantum number (n — 0, 1, 2, . . .). 
After substituting a, 7 , C 2 into Eq.(3.21), we obtain 

Aj . 192ct 4 Z(Z+l) 4V„ 2 a 2 -. 
A/ 1 H 53 p72 1 

n' = -n+-¥ ^— — . (3.22) 

From Eq.(3.20), we find 
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Because for the bound states e > 0, we get 

n' > 0. (3.24) 

If n! > 0, there exist bound states, otherwise, there are no bound states at all. By using 
Eq.(3.21) this relation can be recast into the form 

/-, . 1922(Z+l)a 4 4V 2 a 2 -, 

\ ' W Wc? — 

< n < + ^_ ll (3.25) 

i.e. it gives the finite coupling value. 

The condition 7 > gives the define coupling value for the potential depth V 

o<Vo< 4fe«vWTT) (3 26) 

From Eqs. (3.23) and (3.9), we obtain 



/r,-_vb/i i(i+i)(Ci+c 2 ) , , 



2„ 2 1 ^1(1+1) 



Co 



_o 

^2" 



\ ,0 4v n a 



_S2_ ( -1 _ Z(/+l)(Ci+C 2 ) 
4a 2 1 n , 



(3.27) 



Thus, substituting the expressions of a, C , Ci, C 2 and n' into Eq.(3.27), one can find the 
energy eigenvalues E nt 



( 



TP - v o 



\ 



1 - 



32l(l+l)a 3 



R 



1+ ^ W^~ ln ~ X ) + ~h^ 



± 



1 1921Q+1K _ 4^ 2 a 2 _ _ 1 

Kg fi 2 C 2 




m 2 c 2 + fiWlfl_4a + 12^ 
fl 2 V «0 ~g[ 



192i(i + l)a 4 iV^a 2 



U 



■2n-l 



4V 2 a2 



S 2 

16a 2 



/ 



V 



32Z(Z+l)a 3 



1921(l + l) a 4 4V^ 4V-a^ 
S3 ft 2 "?" I 



(3.28) 



If all two conditions (3.25) and (3.26) are satisfied simultaneously, the bound states exist. 
From Eq.(3.25) is seen that if I = 0, then one gets n < 0. Hence, the Klein-Gordon equation 



8 



for the standard Woods-Saxon potential with zero angular momentum has no bound states. 
For larger values of V (V > 4fiac V^( ;+ jj^ ^ e condition (3.25) is not satisfied. Therefore, no 
bound states exist for these values of V . 

According to Eq.(3.28) the energy eigenvalues depend on the depth of the potential Vo, 
the width of the potential Rq, and the surface thickness a. Any energy eigenvalue must be 
less than Vo- If constraints imposed on n and Vo are satisfied, the bound states appear. 
From Eq.(3.26) is seen that the potential depth increases when the parameter a increases, 
but the parameter R is decreasing for given / quantum number and vice versa. Therefore, 
one can say that the bound states exist within this potential. Thus, the energy spectrum 
Eq.(3.28) are limited, i.e. we have only the finite number of energy eigenvalues. 

The binding energy of a bound Klein-Gordon particle is defined as [13]: 

E b = E nl -m c 2 . (3.29) 

By using the empirical values r = 1.285/m and a = 0.65/m taken from Ref.[14] the 
potential depth V = (40.5 + 0.13A)MeV and the radius of the nucleus R = tqA 1 ^ fm are 
calculated for the atomic mass number of target nucleus A = 40; 56; 66; 92; 140; 208 and 
pions with mass m c 2 = 139.570MeV. In the Table 1. energies of the bound states obtained 
numerically for the spherical standard Woods-Saxon potential for some values of I and n are 
given. 

In addition, we have seen that there are some restrictions on the potential parameters in 
order to obtain bound state solutions. We also point out that the exact results obtained for 
the standard Woods-Saxon potential may have some interesting applications for studying 
different quantum mechanical and nuclear scattering problems. 

However, in the non-relativistic limit with the mapping E^p — m c 2 — > E^ R ^ , according 
to Appendix A for bound state energy eigenvalues, we obtain [8]: 



h?l(l+l) ( -i , 12a 2 



2moa 2 



rp _ rri(t-t-L) I ~t _ 



I | l»2i(! + l)a 4 r, n ^ 1 ^ m V a 2 4i(i + l)a 3 ^ 



_|_ in V {) a 2 



(3.30) 



Now, we are going to determine the radial eigenf unctions of this potential. Having sub- 
stituted 7r(z) and cr(z) into Eq.(2.4) and then solving first-order differential equation, one 
can find the finite function <3>(z) in the interval (0, 1) 



9 



$( z ) = z e {l-z)^ e2 - p2W 



(3.31) 



It is easy to find the second part of the wave function from the definition of weight 
function: 

p(z) = z 2 * (1 - zfV' 2 -^ 2 , (3.32) 
and substituting into Rodrigues relation (2.4), we get 



y n (z) = Bnz-* (1 - z)- 2 ^^ ^- z n+2 < (1 - Z y+*J*=P* 



dz n 



(3.33) 



where B n = ^ is the normalization constant [15]. Then, y n is given by the Jacobi polyno- 
mials 

( 2e,2 x /e 2 -/3 2 + 7 2 ) 

y n (z) = P^ V J (l-2z), (3.34) 

where 

P^\l - 2z) = (1 - z)^ z n+a (1 - z) n+f3 \ . 

n v 1 n\ dz n I v 1 J 

The corresponding u n i(z) radial wave functions are found to be 

1 ) = c^(i- z) ^^pt 2Vm ' 



U n l(Z) 



;i-2«), 



(3.35) 



where C„; is the normalization constants determined using f^°[u n i(r)] 2 dr = 1 constraint, i.e. 



r& f 1 2e-l /i \2a/€ 2 -/3 2 +7 2 -1 

aC ni / z (1 - 2) V 

J o 



PI 



2e,2 v /e 2 -/3 2 + 7 i 



;i - 2*) 



tfc = 1 



IV. CONCLUSION 



In this paper, we have analytically calculated energy eigenvalues of the bound states 
and corresponding eigenfunctions in the new exactly solvable Woods-Saxon potential. The 
energy eigenvalue expression for Woods-Saxon potentials is given by Eq.(3.28). As it should 
be expected, for any given set of parameters Vo,i?o and a, the energy levels of standard 
Woods-Saxon potential are positive. The obtained results are interesting for both theoretical 
and experimental physicists, because they provide exact expression for energy eigenvalues 
and corresponding eigenfunctions. 



10 



Appendix A 



After some simple transformations Eq.(3.27) takes the form 

p(«) _ V ( i l(l+l)(C 1 +C 2 ) \ , 9 1 -, , fi 2 t(i+l)C 

4/L 2 ^+l)C \/ 2 a 2 ft 2 n' 2 f 1 _ Z(m)(Ci+C 2 ) 

fc 2 4 2rf7 /2 i 4v , 2a2 \ Arnica 2 I ,/ , 4V 2 a2\ 

ft 2 c 4 "*oflo(ra' 2 +-R%2- ) \ « 2 ( v ™' 2 + -p ! ^-J 

Expanding, finally, the above eigenvalue in a series of powers of \ yields 




(Al) 



jp( R ) — _Vo 
^nl — 2 



4l(l+l)(d+C 2 ) 



2h?c*( Jl+ 19 W> a4 -2n-l 



4Z(m)(Ci+c 2 ) 

« 2 (y^ 2 raz_ 2n . 



i m r 2 fi , gjO+ljOg 

^(^TiK-2n-l) 2 



/i 2 

2 32m 2 c 



= -M> 4- 

2 



2l(l+l)(C 1 +C 2 )V 



a 2 . 1 



192i(i + l)a 4 



-2n-l 



2 , fi 2 I(i+l)C (1 
m ° C + 2mofig 

ft 2 f(m)(Ci+C 2 ) 
4moa 2 a 2 



1 h 2 
16 2moa 2 



2m V , 2 a 2 



n 2 



fi 2 



32m 



fi 2 f 2 (t+i) 2 (Ci+c 2 ) 2 _ 2 , fi 2 f(t+i)c vb 

/ , r"2" -r --- — ITLqC -\ p2 7T 



2m a 2 a 4 



1 + 



(^l+ 192 ' ( g 1)a4 -2n-l 

) 



192Z(Z+l)a 4 



2n- 1 



2m R 2 1 

/ 2m y- a 2 i^+lXCj+C^ y 

n 2 V P ^ y 

2moa 2 



(A2) 



fi 2 z(/+i)(Ci+c 2 ) , _ 2 , g 2 j(m)Co _ ft 2 

4m a 2 a 2 ~~ m c ~i 2m R 2 2m a 2 



l+ 192 '<j+ 1)a4 -2n-l 



^ ( ^ 1 + 12 Hfi«l 2 l_ 2n _ 1 



1 / 2m Voa 2 _ KWQ(Ci+C^ \ , V 

2 ^ J? a 1 J l 



( 2m V a 2 i(i+l)(C 1 +C 2 \' 

— T? I 



l+ 192 '"+ 1)a4 -2n 



-)■ 



ft 2 *(/+l)Co _ fi 2 
2mo-R ( 2 2moa 2 



H = m c 2 + 

n 2 



1 ( U , 192,(1+1)* 2n ^ 3^!^ 



Denoting, as usual, the energy eigenvalues in the non-relativistic case by — m c 2 



E^ R \ we have 



(JVJi) _ ft 2 z(m)Co 



2m -Rj" 



?l 2 



2moa 2 



n 2 



1 + 



192l(Z+l)a 4 



- 2n - l) 



+ 



2mQ Vq a 

' h 2 



^ 



19 2 I(I + l)a 4 



(A3) 



2n-l 
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Thus, substituting the expressions of Co, C\ and C2 into the above energy eigenvalues 
we obtain 



E. 



(NR) 
nl 



h 2 

2m R 2 



I 4a 

Ro 



12a 2 



a 2 



192?(/+l)q 4 



- 2n - 1 + 



S3 - 



l+ 192 '"+ 1 )° 4 -2n-l 



or 



(l + 



1+ 



192i(i + l)a 

ST" 



-2n-l 



2m i? 2 , 
2 



12a 2 \ 



■i V|)a- 4i(i+l)a 

" 2 



2moa 2 



16 



ll)2?(i + l)a 4 



-2n-l 



+ 



mpViia 
ft 2 
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A 


Rq, fm 


V ,MeV 


n 


I 


E b ,MeV 


40 


4.3946 


45.70 





1 


-107.8777 


56 


4.9162 


47.78 





1 


-127.5238 


56 


4.9162 


47.78 





2 


-17.5985 


66 


5.1930 


49.08 





2 


-50.3359 


92 


5.8010 


52.46 





2 


-101.8967 


140 


6.6724 


58.70 





3 


-92.5327 


208 


7.6136 


67.54 





4 


-105.0865 


208 


7.6136 


67.54 





5 


-33.6014 



TABLE I: Energies of the bound states for the Woods-Saxon potential for different values of n, I 
calculated using Eqs.(3.25), (3.26), (3.28) and (3.29). 
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corresponding eigenfunctions are obtained on the various values of the quantum numbers n and I. 
The non-relativistic limit of the bound state energy spectrum was also found. 
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I. INTRODUCTION 



An analytical solution of the radial part of the Klein-Gordon equation is of high im- 
portance spinless in relativistic quantum mechanics, because the wave function contains all 
necessary information for full description of a quantum system. There are only few poten- 
tials for which the radial part of the Klein-Gordon equation can be solved explicitly for all 
n and I. So far, many methods were developed, such as supersymmetry (SUSY) [1,2] and 
the Pekeris approximation [3-8], to solve radial part of the Klein-Gordon equation exactly 
or quasi-exactly for I ^ within these potentials. 

The one- dimensional Klein-Gordon (KG) equation is investigated for the PT -symmetric 
generalized Woods-Saxon (WS) potential [9] and Hulthen [10] and is solved by using the 
Nikiforov-Uvarov (NU) method which is based on solving the second-order linear differential 
equations by reduction to a generalized equation of hypergeometric type. 

The radial part of the Klein-Gordon equation for the Woods-Saxon potential [11] cannot 
be solved exactly for I ^ 0. But in Refs.[9, 12-14] authors as in Refs.[15, 16] because of errors 
made by them in application of the NU method to investigate the Woods-Saxon potential 
obtained wrong results. It is well known that the Woods-Saxon potential is one of the 
important short-range potentials in physics. And this potential was applied to numerous 
problems, in nuclear and particle physics, atomic physics, condensed matter and chemical 
physics. Therefore, it would be interesting and important to solve the radial of the Klein - 
Gordon equation for Woods-Saxon potential for I ^ 0, since it has been extensively used to 
describe the bound and continuum states of the interacting systems. 

In this work, we solve the radial part of the Klein-Gordon equation for the standard 
Woods-Saxon potential using NU method [17], and obtain the energy eigenvalues and cor- 
responding eigenfunctions for any I states. 

II. NIKIFOROV-UVAROV METHOD 

The Nikiforov-Uvarov (NU) method is based on the solutions of general second-order 
linear equations with special orthogonal functions. It has been extensively used to solve 
the nonrelativistic Schrodinger equation and other Schrodinger-like equations. The one- 
dimensional Schrodinger equation or similar second-order differential equations can be writ- 
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ten using NU method in the following form: 

V{z) + ^\z) + ^(z) = 0, (2.1) 

where a(z) and a(z) are polynomials, at most second-degree, and t(z) is a first-degree 
polynomial. 

Using Eq(2.1) the transformation 

,/>(z) = Hz)y(z) (2.2) 
one reduces it to the hypergeometric-type equation 

a(z)y" + r(z)y' + Xy = 0. (2.3) 
The function <&(z) is defined as the logarithmic derivative [10] 

vu) «,) (24) 



$(z) <T{zy 

where tt(z) is at most the first-degree polynomial. 

The another part of ip(z), namely y(z), is the hypergeometric-type function, that for fixed 
n is given by the Rodriguez relation: 

D Jn 

VnW = rfdd^WzMz)}, (2-5) 

where B n is the normalization constant and the weight function p(z) must satisfy the con- 
dition [12] 

£(a(z)p(z))=r(z)p(z), (2.6) 

with t(z) = t(z) + 2n(z). 

For accomplishment of the conditions imposed on function p(z) of the classical orthogonal 
polynomials, it is necessary, that polynomial t(z) becomes equal to zero in some point of an 
interval (a, b) and derivative of this polynomial for this interval at a(z) > will be negative, 
i.e., t'(z) < 0. 

The function tt(z) and the parameter A required for this method are defined as follows: 



/ ~ / / / ~\ 2 

O — T , I <T — T 



*(*) = — J + ^ ( 2 - 7 ) 
\ = k + n'(z). (2.8) 



On the other hand, in order to find the value of k, the expression under the square root 
must be the square of a polynomial. This is possible only if its discriminant is zero. Thus, 
the new eigenvalue equation for the Eq.(2.1) is [17]: 

A = A n = -nr> - ^_l2 a ", ( n = 0, 1, 2, ...), (2.9) 

After the comparison of Eq.(2.8) with Eq.(2.9), we obtain the energy eigenvalues. 

III. SOLUTIONS OF THE KLEIN-GORDON EQUATION WITH THE WOODS- 
SAXON POTENTIAL 

The standard Woods-Saxon potential [11] is defined by 

V(r) = - V ° (r - RoV a«flo, 0<r<oo. (3.1) 
1 + exp ( z -^f a ) 

This potential was used for description of interaction of a neutron with a heavy nucleus. 
The parameter Rq is interpreted as radius of a nucleus, the parameter a characterizes thick- 
ness of the superficial layer inside, which the potential falls from value V = outside of a 
nucleus up to value V — — V inside of a nucleus. At o = 0, one gets the simple potential 
well with jump of potential on the surface of a nucleus. 

In the spherical coordinates, the stationary Klein-Gordon equation with Woods-Saxon 
potential is [18] 



-h 2 c 2 \ -- (r 2 - ] + 
r 2 dr \ dr 



i d /a 

r 2 sin 9 89 \ 89 



8 2 



r 2 sin 9 dtp 2 



iP=[(E-V(r)) 2 -m%c 4 } ^. 

(3.2) 



where mo is the rest mass of a scalar particle. 

The terms in the square brackets with the overall minus sign are the dimensionless angular 
momentum squared operator, L 2 . Defining ip(r,9,(p) = R(r)Y(9,(p), we obtain the radial 
part of the Klein-Gordon equation [18] 



d 2 R{r) 2 dR(r 



+ — 



+ 



[E - Vy - m 2 c 4 1(1 + 1 



h 2 c 2 



dr 2 ' r dr 
where / is the angular momentum quantum number. 
After introducing the new function: 



R(r) = 0. 



(3.3) 



u(r) = rR(r), 



Eq.(3.3) takes the form: 



d 2 u(r) 
dr 2 



(E - V) 2 - m 2 c 4 1(1 + 1) 



Vc 2 r 2 U ^ = °- (3 ' 4) 

The Eq.(3.4) has the same form as the equation for a particle in one dimension, except 
for two important differences. First, there is a repulsive effective potential proportional to 
the eigenvalue of h 2 l(l + l). Second, the radial function must satisfy the boundary conditions 
u(0) = and u(oo) = 0. 

If in Eq.(3.4) introduce the notations: 

_ r - R _ Ro 

x , ct , 

Ro a 

then the Woods-Saxon potential is given by the expression: 

Vws = ~r— ^7 r- 

1 + exp(o;xJ 

It is known that the radial part Klein-Gordon equation cannot be solved exactly for this 
potential at the value / ^ using the standard methods as SUSY and NU. From Eq.(3.4) 
it is seen that the centrifugal potential Vi(r) = is inverse square potentials and as 

a result, this equation (3.4) cannot be solved analytically. Therefore, in order to solve 
this problem we can take the most widely used and convenient for our purposes Pekeris 
approximation. This approximation is based on the expansion of the centrifugal barrier in 
a series of exponentials depending on the internuclear distance, taking into account terms 
up to second order, so that the effective / dependent potential preserves the original form. 
It should be pointed out, however, that this approximation is valid only for low vibrational 
energy states. By changing the coordinates r = Rq(1 + x), the centrifugal 

potential is expanded in the Taylor series around the point x = (r = R ) 

Tr/ , h 2 l(l + l) h 2 l(l + l) 1 r . n n 2 A o v 

VAr) = — i = — v ' -. = 5 (l - 2x + 3x 2 - Ax 3 + . . .) , 3.5 

v ; 2m r 2 2m R 2 (1 + x) 2 v J> \ J 

where 5 = 

According to the Pekeris approximation, we shall replace potential V/(r) with expression 
[7,8]: 

V7(r) = 5 ( C + + , (3.6) 

V 1 + expax (1 + expaxf) 
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where Co, C\ and C 2 are some constants. 

In order to define the constants C , C\ and C 2 , we also expand this potential in the Taylor 
series around the point x = (r = R ): 



V?(x) = 5 



C + — + — ((7i + C 2 ) x H C 2 x 2 H (Ci + C 2 ) C 2 x 4 + 

2 4 / 4 v ; 16 48 v ; 96 



(3.7) 

Comparing equal powers of x Eqs.(3.5) and (3.7), we obtain the constants Co, C\ and C 2 
[7,8]: 

4 12 8 48 48 

(-^0 — 1 1 o J W — 9 J <^2 — n ■ 



a cr 



a cr 



cr 



Instead of solving the radial part Klein-Gordon equation for the centrifugal potential 
Vi(r) given by Eq.(3.5), we now solve the radial part Klein-Gordon equation for the new 
centrifugal potential V^*(r) given by Eq.(3.6) obtained using the Pekeris approximation. 
Having inserted this new centrifugal potential into Eq.(3.4), we obtain: 



dr 2 



+ 



E't-m^C 4 1(1 + 1) C } \h 2 c2 «2 a 2 J y h 2 c 2 a 2 a 2 



h 2 c 2 



at 2 a? 



1 + e~ 



(l + e^)- 



u = 0. (3.8) 



We use the following dimensionless notations: 



e 2 = 



(E 2 - m 2 c 4 ) a 2 I (I + 1) C \ _ 2 _ 2EV a 2 l(l + l)d , 

; P 



h 2 c 2 



cr 



h 2 c 2 



cr 



5 7=- 



V 2 a 2 l(l + l)C 2 



h 2 c 2 



(3.9) 



cr 



with real e > (E 2 < m^c 2 ) for bound states; /3 is real, 7 is real and positive. 
If we rewrite Eq.(3.8) by using a new variable of the form [4]: 

r-.Ro 



we obtain: 



z — I 1 + exp 



... . 1 — 2z ,. . — e 2 + B 2 z — 7 2 z 2 
{Z) + z(l-z) {Z) + 7Z7T I\v! u ( z ) = °> (° < z < !)> 



(z(l-z)) s 



(3.10) 



with r(z) = 1 — 2^;; cr(z) = z(l — z); a(z) = — 
In the NU- method the new function n(z) is: 



2^2 



P 2 z - -f 2 z 



ir(z) = ±^e 2 + (k- (3 2 )z-(k- 1 2 )z 2 . 



(3.11) 



The constant parameter k can be found employing the condition that the expression 
under the square root has a double zero, i.e., its discriminant is equal to zero. So, there are 
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two possible functions for each k: 

j (e - y/e 2 - /3 2 + 7 2 ) z - e, for k = f3 2 - 2e 2 + 2e^e 2 - (3 2 + 7 2 , 
^ ± I ( e + " /? 2 + 7 2 ) ^ " e, /or jfc = /? 2 - 2e 2 - 2 ev /e 2 - f3 2 + 7 2 . 

According to the NU-method, from the four possible forms of the polynomial Tt(z) we 
select the one for which the function t(z) has the negative derivative and root lies in the 
interval (0,1). Therefore, the appropriate functions tt(z) and t(z) have the following forms: 

n(z) = e - (e + ^e 2 - /3 2 + 7 2 ) z, (3.13) 



t(z) = 1 + 2e - 2 (l + e + ^e 2 - (3 2 + 7 2 ) z, (3.14) 

and 



k = p 2 - 2e 2 - 2e^e 2 - P 2 + 7 2 . (3.15) 
Then, the constant A = k + n'(z) is written as: 

A = f3 2 - 2e 2 - 2 ev /e 2 - /3 2 + 7 2 - e - ^/e 2 ~ /3 2 + 7 2 - (3.16) 

An alternative definition of X n (Eq.(2.9)) is: 

A = A n = 2 (e + v 7 ^ - /? 2 + 7 2 ) n + n(n + 1). (3.17) 

Having compared Eqs.(3.16) and (3.17): 

/ 3 2 _2e 2 -2e v /e 2 - /3 2 + 7 2 -e- ^/e 2 - /3 2 + 7 2 = 2 (e + ^ t 2 - (5 2 + 7 2 ) ra+ra(n+l), (3.18) 

we obtain 

e + ^e 2 - /3 2 + 7 2 + n + \ - ^±±^t = , (3.19) 



or 



e+ v / e 2 -/3 2 + 7 2 -n' = 0. (3.20) 

Here 

(3.21) 

n being the radial quantum number (n — 0, 1, 2, . . .). 
After substituting a, 7 , C 2 into Eq.(3.21), we obtain: 



/ 1 i 192a 4 l(l+l) 4V 2 a 2 _ -. 



From Eq.(3.20), we find: 



e = I L' + ^— ^ ) . (3.23) 



Because for the bound states e > 0, we get: 

n > 0. (3.24) 

This physically interesting point imposes that n' > 0, which certainly determines the 
number of physically meaningful bound states for a deep potential appearing near the sur- 
face. If n! > 0, there exist bound states, otherwise, there are no bound states at all. By 
using Eq.(3.22) this relation can be recast into the form: 

L ~ 1922(Z+l)a 4 _ 4V 2 a2 _ ^ 

< n < ± R " — (3.25) 

i.e. it gives the finite coupling value. 

The condition 7 > gives the define coupling value for the potential depth Vo: 

0< V < (3.26) 

Rq 

From Eqs. (3.23) and (3.9), we obtain: 



Enl = I 1 - W+W]^ I ± 



2 a 2 



cn 



,2„ 2 +ft 2 'q+i), 



\ 12 , iV Q a ' 



4a 2 



1 - 



z(m)(Ci+c 2 ) 



a 2 (n' 2 +- 



(3.27) 



Thus, substituting the expressions of a, C , Ci, C 2 and n' into Eq.(3.27), one can find the 
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energy eigenvalues E n i 



Era = - 



V 



\ 



1 - 



32l(l+l)a 3 



\ 



1+^-^-2" L| 



1 + m ^_4^_ 2n _ 1 



R 



h 2 

16a 2 




± 



:(f+l) 4a. , 12a 2 \ 



192i(i + l)a 4 4V ( fa 2 



h' 2 c' 2 



"2n-l + 1 £t- 



(3.28) 



V 



32Z(7+l)a 3 



1+ 



If all two conditions (3.25) and (3.26) are satisfied simultaneously, the bound states exist. 
From Eq.(3.25) is seen that if I = 0, then one gets n < 0. Hence, the Klein-Gordon equation 
for the standard Woods-Saxon potential with zero angular momentum has no bound states. 
For larger values of V (V > ihac \ / ^( l + }l^ ^ e condition (3.25) is not satisfied. Therefore, no 
bound states exist for these values of V . 

According to Eq.(3.28) the energy eigenvalues depend on the depth of the potential Vq, 
the width of the potential Rq, and the surface thickness a. Any energy eigenvalue must be 
less than Vq. If constraints imposed on n and Vq are satisfied, the bound states appear. 
From Eq.(3.26) is seen that the potential depth increases when the parameter a increases, 
but the parameter R is decreasing for given / quantum number and vice versa. Therefore, 
one can say that the bound states exist within this potential. Thus, the energy spectrum 
Eq.(3.28) are limited, i.e., we have only the finite number of energy eigenvalues. 

The binding energy of a bound Klein-Gordon particle is defined as [18]: 



E b = E nl - m c 2 . 



(3.29) 



By using the empirical values r = 1.285/m and a = 0.65/m taken from Ref.[19] the 
potential depth Vq = (40.5 + 0.13A)MeV and the radius of the nucleus Ro = roA 1 / 3 /???, are 
calculated for the atomic mass number of target nucleus A = 40; 56; 66; 92; 140; 208 and 
pions with mass m c 2 = 139.570Me1/. In Table 1, energies of the bound states obtained 
numerically for the spherical standard Woods-Saxon potential for some values of I and n 
are given. From Table 1 is seen that for fixed n the energy of the bound states increases 
with increase of I. This means that due to the centrifugal potential Vj(r) in the system the 
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repulsive forces appear. Therefore, in order to compensate this potential the energy of the 
bound state must increase [20]. 

In addition, we have seen that there are some restrictions on the potential parameters in 
order to obtain bound state solutions. We also point out that the exact results obtained for 
the standard Woods-Saxon potential may have some interesting applications for studying 
different quantum mechanical and nuclear scattering problems. 

However, in the nonrelativistic limit with the mapping E^p — m c 2 — > E^ R \ according 
to Appendix A for bound state energy eigenvalues, we obtain [8]: 



r, 2 



2m^a 2 



p _ h 2 l(l+l) ( , 12a 2 



m t )Vi)(i 



(3.30) 



Let us note, that in work [12] the Klein-Gordon and Schrodinger equations with Woods- 
Saxon potential for I ^ states was solved, the bound state energy eigenvalues and the 
corresponding eigenfunctions were found. But in the work [12] errors were made in applica- 
tion of the NU method to Schrodinger equation, which led, as in papers [15, 16] to wrong 
predictions for the bound energy eigenvalues and the corresponding eigenfunctions. Such 
errors were made also in the works [9, 13, 14]. 

Now, we are going to determine the radial eigenfunctions of this potential. Having sub- 
stituted ir(z) and a{z) into Eq.(2.4) and then solving first-order differential equation, one 
can find the finite function <3>(z) in the interval (0, 1) 

$(z) = z* {1 - z)^ 2 -^ 2 . (3.31) 

It is easy to find the second part of the wave function from the definition of weight 
function: 

p(z) = z 2 * (1 - z)V< 2 -0 2 +7 2 5 (3.32) 

and substituting into Rodrigues relation (2.5), we get 

d 



n (z) = B n z- 2 * (l - z yW^-^ 1_ 3*+* (1 _ z)n+ 2^qvw 



(3.33) 



where B n = ^ is the normalization constant [21]. Then, y n is given by the Jacobi polyno- 
mials 

(2e,2We 2 -/3 2 +~, 2 ) 

y n (z) = Pr\ V \l-2z), (3.34) 
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where 



1 An r , 

pW\l - 2z) = ^z- a (1 - z)- p ^- z n+a (1 - z) n+f} \ 
n v ' n\ v ' dz n L v 7 J 



The corresponding u n i(z) radial wave functions are found to be 



(z) = C nl z*{l - z)^ 2 -^ 2 F n (: 




(3.35) 



where C n z is the normalization constants determined using j^°[u n i(r)} 2 dr 



1 constraint, 



i.e., 



* I z 2 ^ 1 (1 - z) W*-W-i \ Pj 

J o 



,0 



n 2 



(1-2*) dz = l. 



IV. CONCLUSION 

In this paper, we have analytically calculated energy eigenvalues of the bound states 
and corresponding eigenfunctions in the new exactly solvable Woods-Saxon potential. The 
energy eigenvalue expression for Woods-Saxon potentials is given by Eq.(3.28). The nonrel- 
ativistic limit (3.30) of the bound state energy spectrum was also obtained . As it should be 
expected, for any given set of parameters Vo, Rq and a, the energy levels of standard Woods- 
Saxon potential are positive. The obtained results are interesting for both theoretical and 
experimental physicists, because they provide exact expression for energy eigenvalues and 
corresponding eigenfunctions. 
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Appendix A 



After some simple transformations Eq.(3.27) takes the form: 

p(«) _ V ( i l(l+l)(C 1 +C 2 ) \ , 9 1 -, , fi 2 t(i+l)C 

4/L 2 ^+l)C \/ 2 a 2 ft 2 n' 2 f 1 _ Z(m)(Ci+C 2 ) 

fc 2 4 2rf7 /2 i 4V , 2a2 \ Arnica 2 I ,/ , 4V 2 a2\ 

ft 2 c 4 "*oflo(ra' 2 +-R%2- ) \ « 2 ( v ™' 2 + -p ! ^-J 

Expanding, finally, the above eigenvalue in a series of powers of \ yields: 




(Al) 



^nl — 2 



Al{l+\){Ci+C 2 ) 



+ ™ r 2 fi , fi 2 f(|+l)Co 



S 2 



2fi 2 c 2 (^/l+ 192i( ^ 1)a4 "-2n-l^) 

j 4Z(m)(Ci+C 2 ) 

Q2 (\/ 1+ 192 ' ( H4 1)a4 : - 2 ^- 1 ) 2 
2m V 2 a 2 



2 32m 2 c 2 a 2 



1 + 



192l(l+l)a 4 
Ft 4 



2n - 1 x 



)' 



-= 4- 
2 



2Z(m)(c 1 +c 2 )Vb 



a 2 . 1 



W 2n-l 



2 , h 2 «l+l)C 
h 2 l(l+i){Ci+C 2 ) 

Amoa 2 a 2 



1 h 2 



n 2 



32m 



16 2moa 2 



fi 2 / 2 q+i) 2 (Ci+C2) 2 



I f 192 ^ 1)a4 - 2n - 1 ) f 



+ ~ m 0C + 2mo ^ ~ T (A2) 



1 + 



192;(Z+l)a 4 _ 

^ *' b 1 ) 2m a 2 



( 2m Q V a 2 i('+i)(Ci+C 2 ) 



In- 1 



s 2 V ffr 



/> 2 z(m)(Ci+c 2 ) , _ _ .2 , ft 2 ^+i)Co _ h 2 

4m a 2 a 2 ^ ^ 2m a R 2 2m a 2 



( ] j 1+ljm i^ I - 2n - 1 ) 



1 / gmpjW _ Z(f+l)(Ci+C 2 ) A , 
2^ 7? a 2 y i " 



/ 2m V a 2 l(l + l)(C 1 +C 2 ) ^ " 

I — T? ^ 



( 



l+ 192i(i + 1)a4 -2n 



h 2 l(l+l)C _ n 2 

2mo-R 2 , 2moa 2 



|(^/ 1+ ~l_ 2l ,_ 1 ) + ^ 



H = m c z + 

"I 2 



2m V a* i(i + l)(C 1 +C 2 ) 
2 ^ 



l+ 192i( '+ 1)a4 -2n-l 
R 



Denoting, as usual, the energy eigenvalues in the nonrelativistic case by E^p — m c 2 — >■ 
E^ R \ we have: 



(NR) _ h 2 l(l+l)C 



nl 



2m Rj~ 



h 2 



2m(ja 2 



1 + ^^-27^-1 



+ 



2mQ Vga 

' h 2 



l(l+l)(C 1 +C 2 ) 

TV 



19 2 I(I + l)a 4 
«0 



(A3) 



2n-l 
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Thus, substituting the expressions of Co, C\ and C2 into the above energy eigenvalues 
we obtain: 



E. 



(NR) 
nl 



h 2 

2m R 2 



I 4a 

Ro 



12a 2 



a 2 



192?(/+l)q 4 



- 2n - 1 + 



S3 - 



l+ 192 '"+ 1 )° 4 -2n-l 



or 



(l + 



1+ 



192i(i + l)a 

ST" 



-2n-l 



2m i? 2 , 
2 



12a 2 \ 



■i V|)a- 4i(i+l)a 

" 2 



2moa 2 



16 



ll)2?(i + l)a 4 



-2n-l 



+ 



mpViia 
ft 2 
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A 


Rq, fm 


V ,MeV 


n 


I 


E b ,MeV 


40 


4.3946 


45.70 





1 


-107.8777 


56 


4.9162 


47.78 





1 


-127.5238 


56 


4.9162 


47.78 





2 


-17.5985 


66 


5.1930 


49.08 





2 


-50.3359 


92 


5.8010 


52.46 





2 


-101.8967 


140 


6.6724 


58.70 





3 


-92.5327 


208 


7.6136 


67.54 





4 


-105.0865 


208 


7.6136 


67.54 





5 


-33.6014 



TABLE I: Energies of the bound states for the Woods-Saxon potential for different values of n, I 
calculated using Eqs.(3.25), (3.26), (3.28) and (3.29). 
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